High-accuracy approximation of binary-state dynamics on networks 
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Binary-state dynamics (such as the susceptible- infected-susceptible (SIS) model of disease spread, 
or Glauber spin dynamics) on random networks are accurately approximated using master equations. 
Standard mean-field and pairwise theories are shown to result from seeking approximate solutions of 
the master equations. Applications to the calculation of SIS epidemic thresholds and critical points 
of non-equilibrium spin models are also demonstrated. 

PACS numbers: 89.75.Hc, 64.60.aq, 89.75.Fb, 05.45.-a 

Dynamical processes running on complex networks are used to model a wide variety of phenomena [l], 0] ■ Examples 
include spreading of diseases or opinions through a population 0, 0| , neural activity in the brain [j| , and cascading 
bank defaults in a financial system Q. The structure of the underlying network (e.g., its degree distribution) may 
strongly influence the dynamics and determine critical values of parameters (e.g., the critical temperature of the Ising 
spin model [?] , or the epidemic threshold for disease-spread models 0, H| ) . Accurate prediction of dynamics and 
critical points on networks of arbitrary degree distribution thus remains an important unsolved problem [l|. 

Mean-field theories (MF) are relatively simple to derive and can be quite accurate for dynamics on well-connected 
networks (Io| . However, on sparse networks, or close to critical points, MF theories perform poorly (see, for example, 
Fig. 1 below). Pair- wise approximations (PA), which take into account the states of both nodes at the ends of a 
network edge, improve on MF, but have been derived for fewer dynamical processes (examples are [Tl, 12|). In this 
Letter we demonstrate a tractable master equation approach for binary-state dynamics, with accuracy exceeding both 
MF and PA. We show that PA and MF theories may be derived by seeking approximate solutions of the master 
equations. We write down the explicit PA equations for the general case, thus giving the first derivation of pair- 
wise approximations for a range of dynamical processes. Finally, we use the master equations to calculate critical 
points such as the epidemic threshold for the susceptible-infected-susceptible (SIS) model (or contact process), and 
the critical noise level in the majority- vote model |16j |. 

We consider binary-state dynamics on static, undirected, connected networks in the limit of infinite network size. For 
convenience, we call the two possible states of a node susceptible and infected, as is common in disease-spread models. 
However, this approach also applies to other binary-state dynamics, such as spin models [17|, where each node may 
be in the +1 (spin-up=infected) or the —1 (spin-down=susceptible) state. The networks have degree distribution 
and are generated by the configuration model [2]. Dynamics are stochastic, and are defined by infection and recovery 
probabilities which depend on the degree fc of a node, and on the current number m of infected neighbors of the node. 
Thus Fk^m dt is defined as the probability that a fc-degree node that is susceptible at time t, with m infected neighbors, 
changes its state to infected by time t + dt, where dt in an infinitcsimally small time interval. Similarly, Rk, m dt is 
the probability that a fc-degree infected node with m infected neighbor moves to the susceptible state within a time 
dt. These general infection and recovery probabilities can describe many dynamical processes of interest, see Table U 
for some examples. 
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SIS [13] 
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Voter model [141 


m/fc 


1 Fk,m 


Glauber dynamics [15] 


[l + exp(^-(fc-2m))] _1 


1 -Ffc,m 




f Q ifm<fc/2 




Majority- vote [161 


< 1/2 if m = ft/2 


1 - F k , m 




[l-Q if m > fc/2 





TABLE I: Infection and recovery rates for some examples of binary-state dynamics on networks: k is the node's degree, m is its 
number of infected neighbors. Parameters A and \i are SIS transmission and recovery rates; T and J are the temperature and 
interaction strength for the Ising model; Q is the majority-vote noise parameter. Note T = Glauber dynamics are identical 
to those of the Q — majority-vote model. 
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Approximate master equations for dynamics of this type can be derived by generalizing the approach used in [18 1 
for SIS dynamics, see Appendix A. Let Sk, m (t) (resp. ik,m{t)) be the fraction of fc-degree nodes that are susceptible 
(resp. infected) at time t, and have m infected neighbors. Then the fraction Pk(t) of fc-degree nodes that are infected 
at time t is given by pk(t) — X)m=o * fe > m ~ ^ ~ El=o s k,m, and the fraction of infected nodes in the whole network is 
found by summing over all fc-classes: p(t) — (pk{t)) = J2k Fk Pk{t)- 

The master equations for the evolution of s k>m (t) and ik,m(t) are (see Appendix A): 

4: s k,m = -F k , m Sk,m + Rk, m ik, m ~ P s (k - m)s km + /3 s (k - m+ l)sfc, m _i - 7 s ms fe m + 7 s (m+ l)s fc m+1 , (1) 
at 

jl^.m = -Rk, m ik.rn + F k , m Sk, m ~ (k - m)i k , m + /3 l (k - m + l)i k , m -i - 7 l mi fcim + 7 l (m + l)i fcjm+ i, (2) 

for each m in the range 0, . . . , k, and for each fc-class in the network. The first two terms on the right hand side 
of each equation represent transitions due to infection or recovery of a fc-degree node. The remaining four terms 
account for infection or recovery of a neighbor. The rates (3 s , 7 s , f3 l , and 7* are approximated by tracking the number 
of edges of each type. To calculate (3 s , for example, we count the number of S-S edges (i.e., edges between two 
susceptible nodes) in the network at time t, and then count the number of edges which switch from being S-S edges 
to S-I edges in the time interval dt; the probability /3 s dt is given by taking the ratio of the latter to the former, 

* = ^asaer 1 - ~ - - f - 'is^-sa-' . p = ™ d 

7* = ^£pf m Rk,m lk ( n ^ , see Appendix A for details. 

The master equations (p} and <j2j) , with the time-dependent rates /3 s , 7 s , f3 l and 7* (defined as nonlinear functions 
of Sk,m and ik.m), form a closed system of deterministic equations which can be solved numerically using standard 
methods. Assuming a randomly-chosen fraction p(0) of nodes are initially infected, the initial conditions are Sfc jTO (0) = 
(1 - p(0)) B ktm (p(0)), i k>m (0) = p(0)B ktm (p(0)), where B h , m (q) denotes the binomial factor (? n )q m (l - qf- m : Note 
that the evolution equations are completely prescribed by the functions F k<m and R ki m, and so this method can be 
applied to any stochastic dynamical process defined by transition rates of this type. For the SIS model, equations ([1]) 
and © were derived in [18| |. but were not analyzed as here. 

Figure 1(a) shows the infected fraction p(t) of nodes in the SIS model run on a 3-regular random graph (i.e., 
a Bethe lattice, with P k = 8 k $)- The master equations ([I])-([2]) clearly give a better approximation to the actual 
stochastic dynamics than standard methods (here, the mean-field theory of [3( and the pair-approximation method of 



111 1 1 Qj j — note these are reproduced by equations ((H) and ^ below) . The steady-state infected fraction is plotted as 
a function of the non-dimensional recovery rate p/X in Fig. 1(b). The master equation solutions give a significantly 
better estimate of the epidemic threshold than the standard approximations: we pursue this further below. Figures 
HJa) and HJb) demonstrate that similar conclusions hold for zero-temperature Glauber dynamics [2(| on networks 
with truncated power-law degree distributions and on 3-regular random graphs. Here the comparison is with the 
mean- field theory of [2l[ (see also ((U) below), and the pair approximation from equation ([3]) below. Figure 2(b) shows 
that our approach captures the fact that T = Glauber dynamics on networks can freeze in disordered states; this 
phenomenon is not captured at all by MF [2lj [ . 

For dynamics on a general network, with non-empty degree classes from k = up to a cutoff fc max , the number of 
differential equations in the system (P)-© is (fc max + l)(fcmax + 2), and so grows with the square of the largest degree. 
In certain no-recovery cases (i.e., Rk, m = 0), such as Watts' threshold model [22j], fc-core size calculations [23[, and 




FIG. 1: (a) Infected fraction p(t) in the SIS disease spread model on 3-regular random graphs, with transmission rate A = 1 
and recovery rate p = 1.4. (b) Steady-state fraction of infected nodes as a function of the non-dimensional recovery rate p/\. 
The arrow marks the epidemic threshold predicted from the linearized master equations (top row of Table HIT a)). 
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bootstrap percolation 24], we can show that an exact solution of the master equations is obtained by solving just 
two differential equations (as given in [25j]). For general dynamics, however, some approximation is necessary if it is 
desirable to reduce the master equations to a lower-dimensional system. One possibility is to consider the parameters 
Pk{t) (resp. qk(t)), defined as the probability that a randomly-chosen neighbor of a susceptible (resp. infected) k- 
degree node is infected at time t. Noting that Pk(t) can be expressed in terms of Sk, m as X)m=o ms k,m/ El=o ksk, m , 
an evolution equation for pk may be derived by multiplying equation (TTJ) by m and summing over m. The right- 
hand-side of the resulting equation contains higher moments of Sfc, m , so a closure approximation is needed to proceed. 
If we make the ansatz that Sf. m and ik, m are proportional to binomial distributions: Sk. m ~ (1 — Pk) Bk, m {Pk), 
ikm ~ Pk Bk,m{<lk)i, we obtain the pair approximation (PA), consisting of the 3/c max + 1 differential equations: 



dt Pk 



-pk 



Rk,m.Bk,m{qk) + (1 — Pk) ^ Fk, m Bk,m(Pk) , 



dt Pk 



dt 



(Ik 



771=0 

k 

£[■ 

m— 



Qk 



Fk.mBk.rn (,Pk ) 
Rk,mBk,m{qk) 



Pk 



1 - Pk 
1- Pk 



Pk 



Rk,mBk,m(qk) 
Fk, m Bk,m(Pk) 



+ P s (l-p k )-J s p k , 

+ W(i~q k )-Yqk, 



(3) 



The rates here are given by inserting the binomial ansatz into the general formulas, so that (3 s , for 
- Pk) J2 m ( k - m ) F k, m B k:m {p k )) I ((1 - Pk)k{\ ~Pk))] initial conditions are p k {0) = p k (0) = q k (0) = 



for each fc-class. 
example, is ((1 - 
P(0). 

A cruder, mean-field (MF), approximation results from replacing both p k and q k with u>: Sk, m ~ (1 — p k ) B k , m (uj), 
ik.m ~ Pk Bk.m{u), where u> = (f Pfc) is the probability that one end of a randomly-chosen edge is infected. Using this 
ansatz in the master equations yields a closed system of fc max + 1 differential equations for the fraction p k of infected 
fc-degree nodes: 



d 

dt pk 



-Pk 



E 

771=0 



Rk, m Bk,m{u) + (1 — Pk) ^2 Fk,mBk,rn{u), 



(4) 



771=0 



with p fc (0) = p(0). 

The PA and MF approximations Q and ((3]) yield increasingly simpler systems of equations for any process that 
can be expressed in terms of infection and recovery rates F^ m and Rk, m - For the SIS model, the PA equations 
([3]) are those of Eames and and Kee ling while the MF equations (HJ) are precisely those of Pastor-Satorras and 
Vespignani For the voter model [l4| . the MF equations (dJ reduce to those in [26|, while the PA equations Q lie 
between those of [11] and [13] in terms of complexity. The MF equations ((4]) for zero-temperature Glauber dynamics 
reproduce the mean- field theory of (2l| (in the limit of infinite network size). For this and related non-equilibrium 
spin models, such as the majority- vote model, steady-state PA equations for the special case of 4-regular graphs (i.e., 
Pk = 5k, 4) are derived in [17] . However, to our knowledge, no PA equations such as ([3]) have been derived for these 
dynamics on networks with arbitrary degree distribution Pk. Note also that a coarser type of PA, using the ansatz 
Sfe,m = (1 — Pk)Bk.m(p), ik.m = PkB k.rn.i o) (i.e., with /c-independent parameters p and q) gives the equations recently 
derived in [28j for SIS, and those in [27J for the voter model. See Appendix B for details of this "homogeneous" PA. 




FIG. 2: Infected fraction p(t) (i.e., fraction of +1 spins) for zero-temperature Glauber dynamics on (a) networks with truncated 
power-law degree distribution: Pk oc k~ 2 ' 5 for 3 < k < 20, and (b) 3-regular random graphs. In each case the initial condition 
is p(0) = 0.4. 
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(a) SIS, z-regular 


(b) Majority- vote, PRG 
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5 


6 


0.240 
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10 


(8.63, 9) 


8.97 


9 


10 
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0.275 


0.277 
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TABLE II: (a) Critical values of fi/X for epidemic spread in the SIS model on z-regular graphs. Lower and upper bounds for 
the critical value of fi/\ for the contact process on a tree (defined as the largest value of /z/A for which the infection survives 
forever with positive probability) are from 12911 . Note that the largest eigenvalue of the adjacency matrix for these networks is 
Ai = z, so the method of Prakash et al. [31j gives the same (inaccurate) prediction for the critical value as MF theory, (b) 
Critical value of the noise parameter Q in the majority-vote model on Poisson (Erdos-Renyi) random graphs of mean degree 
(k) = z. Numerical values are from [3y], other values are determined via stability analysis of equations |T}-((3}. 



We briefly highlight another important application of the master equations: the calculation of the epidemic threshold 
for the SIS disease spread model [H, E|. If the seed fraction of infected nodes p(0) is sufficiently small, an appropriate 
linearization of the master equations (HJ) — (EJ determines whether the infected fraction will grow (to epidemic pro- 
portions), or will decay to zero. This reduces the problem to linear stability analysis, and so to the calculation of 
the largest eigenvalue of a matrix (with dimension of order fc^ ax ). In Table HTJa) we show the critical values of the 
parameter fi/X for SIS dynamics on z-regular random graphs calculated in this way, and compare with the explicit 
values predicted by PA [ill EH and MF Q methods ( i.e., z — 1 and z, respectively). Recently it was argued that SIS 
infection can persist indefinitely in networks containing nodes of sufficiently high degree, due to recurring reinfections 
between hub nodes and their neighbors [9j. The master equation formalism does not capture this effect, because the 
definitions of the rates (/3 s , 7 s , etc) use global counts of edge types, and so wash out structural correlations specific 
to the immediate neighborhood of hub nodes. 

Linear stability analysis may also be applied to spin models with up-down symmetry, which have the property 
Rk,m = 1 — Fk.m — Fk^k-m,, and where the magnetization M(t) (the average of all spins in the network) is given by 
M = 2/9—1. Stability analysis of the (disordered) fixed point with p = 1/2 gives the location of critical points marking 
the transition between disordered and ordered phases. Applying this method to Glauber dynamics reproduces the 
results of for the critical temperature of the Ising model. It also accurately approximates numerically-determined 
critical values for non-equilibrium spin models, such as the critical noise Q c in the majority-vote model, see Table 
11(b). 

In summary, we have derived the master equations (d])-© — first introduced for SIS dynamics in [18[ — for general 
binary-state dynamics on networks, and demonstrated that their accuracy supersedes standard MF and PA methods. 
Mean-field and pairwise theories are derived as approximate solutions of the master equations, and equations ([3]) 
explicitly give pair approximations for any dynamics defined by infection and recovery rates Fk, m and Rk, m - Finally, 
we demonstrated the application of the master equations to calculating epidemic thresholds and critical parameter 
values via linear stability analysis, improving significantly on existing MF and PA estimates. 

We anticipate further applications of the master equation approach to the calculation of critical points in opinion 
models and spi n sy stems, and expect possible extensions to include multiple-state dynamics (such as the SIR disease- 
spread model l32j - [35j |). multiple node types [36j |. discrete-time dynamics [37| . and network models with non-zero 
clustering 38l442j|. 

This work was funded by Science Foundation Ireland awards 06/IN. 1/1366 and MACSI 06/MI/005. Helpful dis- 
cussions with Sergey Melnik, Rick Durrett, and Claudio Castellano, and participants at the SAMSI Dynamics on 
Networks workshop are gratefully acknowledged. 



Appendix A: Derivation of Master equations 

[The material in Appendices A and B appeared in the early ArXiv versions of this paper.] 

We consider binary-state dynamics on static, undirected, connected networks in the limit of infinite network size 
(i.e., N —> 00, where N is the number of nodes in the network). For convenience, we call the two possible node 
states susceptible and infected, as is common in disease-spread models. However, our approach also applies to other 



Ik,m-1 set 



Sk,m set 



I km set 



7 l 



Sfc/m+l set 




•5 




X = 



r = 



FIG. 3: Schematic of transitions to/from the Sj^m and Ik,m sets, as described in equations (|12[) through (|20|) . For each set, the 
central (Ego) node is shown along with some of its neighbors: black nodes are infected, white nodes are susceptible. See also 
Fig. 1 of [3- 



binary-state dynamics, such as spin systems .17] , where each node may be in the +1 (spin- up) or the —1 (spin-down) 
state. The networks have degree distribution Pk and are generated by the configuration model [2,143, 44 [. Dynamics 
are stochastic, and are defined by infection and recovery probabilities which depend on the degree fc of a node, and 
on the current number m of infected neighbors of the node. Thus Fk, m dt is defined as the probability that a k-degree 
node that is susceptible at time t, with m infected neighbors, changes its state to infected by time t + dt, where dt in an 
infinitesimally small time interval. Similarly, Rk, m dt is the probability that a fc-degree infected node with m infected 
neighbor moves to the susceptible state within a time dt. These general infection and recovery probabilities can 
describe many dynamical processes on networks. For example, in the susceptible-infected-susceptible (SIS) model of 
disease spread (the contact process) each susceptible node may be infected at a rate A by each infected neighbor, 
and each infected node recovers at a constant rate /i, so the rates Fk m and Rk,m take the form 



F, 



sis 



A m, 



R 



sis 



(5) 



Zero-temperature Glauber dynamics for a spin system on a network [HI, 2^, 21 1 provide another example. Here, each 
node has a spin of +1 or —1 (which we can identify with the infected or susceptible state, respectively). In each 
infinitesimal time step, one node is selected at random (dt = 1/N) and its spin is set to +1 if the local field (sum of 
its neighbors' spins) is positive, to —1 is the local field is negative, and to ±1 with equal probability if the local field 
is zero. Thus, the new spin matches the majority of its neighbors' spins (with a random choice in case of a tie). The 
infection and recovery rates may therefore be expressed as 



pGlauber 



i£m<k/2, 
| ifm = fc/2, 



oGlaubcr 



1 if m > fc/2, 



1 if m< fc/2, 
| if m= fc/2, 
if m > fc/2. 



(6) 
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The voter model 14j has rates given by 



pvoter _ m nvotcr _ fe m (n\ 

r k,m — ~r> "•fe.m — T > \< ) 



since infection of a fc-degree node, for example, occurs by copying one of m infected neighbors out of fc possible choices. 

We now proceed to derive the master equations for dynamics of this type, closely following the approach used in 
[lij for SIS dynamics. Let Sk, m (resp. Ik, m ) be the set of nodes which are susceptible (resp. infected), have degree 
k, and have m infected neighbors. To quantify the size of these sets, define Sf., m (t) (resp. ik, m (t)) as the fraction of 
fc-degree nodes that are susceptible (resp. infected) at time t, and have m infected neighbors. Then the fraction Pk(t) 
of fc-degree nodes that are infected at time t is given by 

k k 
m— m— 

and the fraction of infected nodes in the whole network is found by summing over all fc-classes: 

p{t)=Y,PkPk{t). (9) 

If a randomly-chosen fraction p(0) of nodes are initially infected, then the initial conditions for Sk iTn and ikm are 
easily seen to be 

s fe ,m(0) = (1 - p(0)) B fc , m (p(0)), i fc , m (0) = p(0)£ fc , m (p(0)), (10) 

/A 

where we introduce the convenient notation Bk m(s) f° r the binomial factor \ q m (l — q) k m . Note that we can 

also calculate the number of edges of various types using this formalism. For example, the number of edges in the 
network which join a susceptible node to an infected node (we call these S-I edges for short) can be expressed in two 
equivalent ways: 

fe k 

N ^2 p k^2 ms k , m or N^2Pk^2(k-m)ik, m - (11) 

k m— k m—0 

The first of these expressions, for example, follows from noting that in a sufficiently large network that there are NPk 
nodes of degree fc, of which a fraction Sk, m are susceptible and have m infected neighbors. Each such node contributes 
m edges to the total number of S-I edges. Similar expressions may also be given for the number of S-S and I-I edges 
in the network. We note that the equivalence of the two expressions in is preserved by the evolution equations 
described below. 

Next, we examine how the size of the Sk, m set changes in time. We write the general expression 

Sk,m(t + dt) = Sk,m{t) - W{Sk,m-*h,m) Sfc.m dt + W(7fc, m ->- Sk,m) ifc.m dt 

— W(Sk,m^Sk,m+l) s k,m dt + W(Sk,rn-l — > Sk,m) Sfe,-m-l dt 

—W(Sk,m—>Sk,m-l) Sk,m dt + W(Sk, m+1 — > Sk,m) Sfc,m+1 dt (12) 

to reflect all the transitions whose rate is linear in dt (all other state-transitions are negligible in the dt — > limit), 
see Fig. [3] Here W(Sk, m ~^ Ik,m) dt, for example, is the probability that a node in the Sk, m set at time t moves to 
the Ik, m set by time t + dt. It is clear from the definitions above that 

W(S ktm -*h,m) = F k , m and W(I k , m ^-Sk, m ) = Rk,m- (13) 

A node moves from the Sk, m -i set to the Sk,m set if it remains susceptible, while one of its susceptible neighbors 
becomes infected. Note this means that an S-S edge ch ang es to an S-I edge. If we suppose that S-S edges change to 
S-I edges at a (time-dependent) rate (3 s , we can write [45| 

W{S Km ^S Km+l )=p s {k-m) and W{S Km _ x ->S fe , m ) = /3 s (k - m + 1), (14) 

since nodes in the Sk, m set have k — m susceptible neighbors, while those in the SV-,m-i set have fc — m + 1 susceptible 
neighbors. To calculate /3 s , we count the number of S-S edges in the network at time t, and then count the number of 
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edges which switch from being S-S edges to S-I edges in the time interval dt; the probability /3 s dt is given by taking 
the ratio of the latter to the former, i.e. 

^8 dt = Efc P k ELo( fe " m ) F k,rn S k , m dt ^ 

A similar approximation is used to define 7 s , the rate at which S-I edges change to S-S edges due to the recovery of 
an infected node: 



Sfc p k Sm=o( fc ~ m)R Km i k , 

Efc P fcEm=o( fc - TO )*fc,m 



(16) 



and we then write 

W(S k , m ^S k , m -i) =l s m and W(S k , m +i ->S fc , m ) = Y{™ + !)■ (17) 
Taking the limit dt — > of equation (fT2")) gives the master equation for the evolution of Sfc, m (i) (see Fig. 

4f. s k,m = -Fh, m Sk,m + Rk,mik,m ~ P s (k - m)sfc, m + /3 s (k - m + l)s k , m -i - 7 s msfe, m + 7 s (m + l)sfc, m+ i, (18) 

where m is in the range 0, . . . , k for each fc-class in the network (and adopting the convention s k ,-i = Sk,k+i = 0). 
Applying identical arguments, mutatis mutandis, to the set Ik, m , we derive the corresponding system of equations for 

^k 7 m (^) • 

4^k,m = -Rk, m ik, m + ft,mSt,m - f?{k - m)i k ^ m + I3' l {k —m + l)ifc, m -i - 7*™fc,m + 7*( m + l)i fc)m+1 , (19) 
for m = 0, . . . , k and for each fc-class in the network, with time-dependent rates f3 z and 7* defined though s kiTn and 

^fc,m as 

Efc P k Em=0 TO S fe,m Efe P fc Em=0 TO 

The master equations ([180 and <JX9j) , with the time-dependent rates /3 s , 7 s , /3* and 7* (defined as nonlinear functions 
of Sfc iTO and i k ,m), form a closed system of deterministic equations which, along with initial conditions (fTUj) . can be 
solved numerically using standard methods [4q |. Note that the evolution equations are completely prescribed by the 
functions F k , m an d Rk,m, and so this method can be applied to any stochastic dynamical process defined by transition 
rates F k , m an d R k , m - For the SIS model, equations (TT8l) and (TT9l) were derived in [3] (see also (HI), with additional 
terms to study adaptive rewiring of the network. 



Appendix B: Homogeneous pair approximation 



The pair approximation J3]) derived in the main text is of the type dubbed "heterogeneous PA" in [12J , because the 
system includes variables p k and q k for each degree class k. A more parsimonious set of equations may be derived 
under the assumptions of "homogeneous PA" , wherein the fc-dependence of edge-based variables p k and q k is neglected. 
As discussed in [12J , the reduction in the number of variables typically comes at the cost of reduced accuracy. 

For homogeneous PA, the parameter p(t) (resp. q(t)) is defined as the probability that a randomly-chosen neighbor 
of a susceptible (resp. infected) node is infected. Noting that p(t) can be expressed in terms of s kym as 



, , , _ Efc P k Em=0 171 Sk ,m _ Efc P k Em=0 171 s k,m y , 

~ Y, k PkEi= ksk, m EkPkHi-Pk) ' 



an evolution equation for p may be derived by multiplying equation (|18p by P k m and summing over m and k. The 
right-hand-side of the resulting equation contains higher moments of s ktTn , so a closure approximation is needed to 
proceed. If, similar to the steps yielding equations ([3]), we here make the ans&tz that s^.m and m are proportional 
to binomial distributions: 

Sfe.m ~ (1 - Pk) -Bfe )OT (p), i k , m ~ p k B kim {q), (22) 
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then an equation for dp/dt may be found in terms of only p k , p, and q. Applying the same ansatz to expression 
gives the algebraic relation 

(l-9)w=p(l-w) (23) 

between p and q, where oj = ^2 k -Pf. can be interpreted as the probability that the node at one end of a randomly- 
chosen edge is infected. After some algebra, we obtain the homogeneous pair approximation, consisting of the fc max + 2 
differential equations: 

4:Pk = ~Pk ^2 R k,rnB k ,m(q) + (1 - Pk) ^ Fk,m B k,m(p) for fc = 0, . . . , fc max , 
m m 

= j-^J 7^ ( 1+P_2 t) ^ ~ P^ Fk < mBk < m< ^ - PkRk, m Bk, m (q)) , (24) 
fe m 

along with the algebraic relation (|23|) . and initial conditions /9fc(0) =p(0) = p(0). 

For the SIS model, the homogeneous PA equations (p4|) are identical to those recently derived by House and Keeling 
[28| . while for the voter model equations (f24"|) are equivalent to those in [27| (in the N — > oo limit). Note that on 
z-regular graphs, the heterogeneous and homogeneous pair approximations are identical. 
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